Solitonic transmission of Bose-Einstein matter waves 
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We consider a continuous atom laser propagating through a wave guide with a constriction. Two 
different types of transmitted stationary flow are possible. The first one coincides, at low incident 
current, with the non-interacting flow. As the incident flux increases, the repulsive interactions de- 
crease the corresponding transmission coefficient. The second type of flow only occurs for sufficiently 
large incident currents and has a solitonic structure. Remarkably, for any chemical potential there 
always exists a value of the incident flux at which the solitonic flow is perfectly transmitted. 
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The transport properties of matter confined to small 
structures display distinct quantum effects qualitatively 
different from those observed at macroscopic scales. 
These are grounded on global phase coherence through- 
out the sample and can be, in many cases, understood 
within single particle pictures, without referring to any 
specific details of the system. As a result, they arise in 
many different fields (electronic systems, atomic physics, 
electromagnetism, acoustics) [2-4]. Some examples are 
(weak and strong) localization, Bloch oscillations and 
conductance quantization. 

Recent experimental developments in the physics of 
Bose-Einstein condensation (BEC) of dilute vapor (in 
particular the microchip guiding technique) open up the 
prospect of studying coherent transport phenomena us- 
ing guided atom lasers [5]. Besides, because of the 
extraordinary control over these systems, they offer a 
unique opportunity to go beyond the single particle be- 
havior, and to study specific effects induced by interac- 
tion. In the present article we focus on a simple situation, 
where a BEC matter wave propagates through a guide 
with a constriction [6]. By an adiabatic approximation, 
the three-dimensional flow is reduced to one dimension, 
where the atoms now feel, due to the constriction, a lon- 
gitudinal step-like potential of height Vq. In the absence 
of interaction, the transmission T does not depend on 
the incident current but only on the beam's energy; T 
is always lower than unity and tends to this limit when 
the energy of the beam is large compared to Vq. In the 
following we consider atoms with a repulsive effective in- 
teraction characterized by a scattering length a sc > 0. 
The most salient features of the flow are all at variance 
with respect to the non-interacting case: (i) the trans- 
mission coefficient depends on the current, (ii) at given 
chemical potential, there exists a maximum transmitted 
current above which no stationary flow exists, (iii) at a 
given current, several distinct stationary solutions with 
different T are possible and (iv) for any chemical po- 
tential larger than Vq, there is a particular value of the 
incident current which induces total transmission. 

Consider a continuous atom laser incident on a con- 
striction of a waveguide. Within the adiabatic approxi- 



mation [7] the transverse motion is restricted to the low- 
est transverse eigen-state. The constriction affects the 
longitudinal motion via an effective step-like potential 
whose magnitude is fixed by the ground state energy of 
the transverse Hamiltonian. For a BEC system, the adi- 
abatic approximation implies that the condensate wave 
function can be cast in the form (see Ref. [8]), 

¥(f,t) =ip(x,t)cf)(r ± ;n;x) , (1) 

where ip(x,t) describes the motion along the axis of the 
laser (the beam is flowing along the positive x direction). 
<fi is the equilibrium wave function (normalized to unity) 
in the transverse (rj_) direction. It depends parametri- 
cally on the longitudinal density n(x,t) = J" d 2 rj_| 1 I r | 2 = 
\4>(x,t)\ 2 . The beam is confined in the transverse di- 
rection by a trapping potential V± (rj_ ; x) , which is x- 
dependent in the region of the constriction. Then, the 
longitudinal wave equation reads [8,9] (in units where 
h = m = 1) 

-\d xx t/j + {v {l (x) + e[n{x, t); a]} V = « d t ip . (2) 

In (2), Vj| (x) represents an effective longitudinal poten- 
tial due to the constriction. If, to be specific, we consider 
a transverse harmonic confinement with pulsation uj±_ (x) , 
then Vj|(x) = w±(x) ~lu±(—oo) (energy is measured with 
respect of the ground state energy of the non interact- 
ing transverse Hamiltonian far before the constriction). 
e(n;x) is a nonlinear term describing the mean field in- 
teraction averaged over a transverse slice of the beam. 
One has e(n; x) = 2u)_\_(x) na sc in the low density regime 
(na sc <C 1), and e(n;x) = 2tu±(x)^/na sc in the high den- 
sity regime (na sc ^> 1) [8,9]. 

Our purpose is to determine the transmission 
of steady state solutions of (2) where tp(x, t) — 
exp{— i[it}A(x) exp{iS(x)}, with A and S real functions. 
The density is n = A 2 and the local velocity is v — dS/ dx. 
From Eq. (2) one obtains (i) flux conservation: n(x)v{x) 
is a constant that we denote Joo , and (ii) a Schrodinger- 
like equation for the amplitude: 

-k^ + {^ + ^y^ + 2^} A= ^ A - (3) 
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To define the scattering problem one needs to study the 
asymptotic behavior of the flow far from the constriction. 
Far upstream V\\(x — ► — oo) = and the non-linear term 
in (2) looses its explicit x dependence, taking the simpler 
form e[n(x,t)]. Thus, in this region, (3) admits a first 
integral of the form [9] : 



dAY 
dx J 



+ W[n(x)] = E d 
with W(n) 



(4) 



J 2 

£ (n)+ M n+-2^ , 
2 n 



where e(n) = J™ e(p)dp and E c \ is an integration con- 
stant. Eq. (4) has a simple interpretation in terms of 
classical dynamics. It expresses the energy conservation 
of a fictitious classical particle with "position" A and 
"time" x, moving in a potential W; E c i being the to- 
tal energy of this particle. Eq. (4) is thus integrable by 
quadrature, and the density profile can be deduced from 
the plot of Fig. 1. Small values of E c i — W{n\) corre- 
spond to small density oscillations, whereas the highest 
acceptable value is E c i = W(n 2 ), corresponding to a gray 
soliton. 
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FIG. 1. W as a function of n. n\ and ni are the zeros 
of dW/dn. At a given E c i, the up-stream density (or equiv- 
alently the "position" of the fictitious particle) oscillates be- 
tween timin and n m ax defined by W(n min ) = W(n max ) = E cl . 

In the far down-stream region, Vj| (x — > +oo) also takes 
a constant value Vq = lu±(+oo) — lu±(—oo) > 0. Hence, 
(3) admits, in this region, a first integral analogous of (4) 
where, due to the change in e(n) (resp. s(n)) takes 
a different form which we denote £o( n ) (resp. £o(n)). 
The new form of W(n) is denoted Wo(n) and the new 
constant of integration is E®f 



(5) 



with W (n) = -eo(n) + (ji - V ) n + . 

It follows from general arguments on the dispersion of 
elementary excitations of Eq. (2) that the physically ac- 
ceptable boundary conditions of Eq. (3) correspond to a 
constant far down-stream density (see [9]). The asymp- 
totic x — > +oo density should thus be equal either to m,o 



(we denote this as "case A"), or to n 2 .o (case B); m.o and 
"2,0 _ being the analogous of ni and n 2 of Fig. 1 - are 
extrema of Wq (n) . They are solutions of 



(6) 



In the non- interacting case the term £o(n) is absent from 
Wo(n) which has only one minimum (nio). Case A is 
therefore the only possible solution in non-interacting 
systems. Case B describes new non-perturbative effects 
related to interaction. It corresponds to an asymptotic 
down-stream density which is part of a gray soliton. 

There exists a maximum value J^ ax of J x above 
which Eq. (6) admits no solution: as Joo is increased 
(keeping \i and Vq fixed), the two extrema of Wo(n) move 
toward each other, until they coalesce and disappear. 
This marks the onset of a time-dependent flow. If, to 
be specific, we consider the case e (n) = go n Va , then 



J" 



-M-v ) 
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l/vo 



(7) 



The scattering process is now well defined. It corresponds 
to the matching between two asymptotic densities de- 
scribed by the classical motion of a particle of energy E c i 
in a potential W(n) at x — > — oo, and of energy E^ in 
a potential Wo(n) at x — > +oo (with E^ either equal to 
Wo(tii,o) or to Wo(n 2 .o))- Eq. (3) being non-linear, an 
important question is how to properly define a transmis- 
sion and a reflection coefficient; i.e., is it possible to dis- 
entangle an incident and a reflected wave in the upstream 
flow ? We follow here an approach closely related to usual 
experimental set-ups, and choose to work with an inci- 
dent and a reflected beam which can be approximated 
by plane waves. This corresponds to a regime where 
Eqs. (3,4) can be linearized in the far upstream region. 
In this regime, for x — > — oo, we write n(x) = m + 5n(x) 
and expand nW(n) to second order in 5n. Then Eq. (4) 
leads to 



dSn 
dx 



+ Kj5n 2 =8{n 1 + Sn)[E c i-W(n 1 )] , (8) 



where k\ — A(v 2 — c\), v\ = Joo/«i being the average ve- 
locity of the up-stream beam and c\ = [ni(de/ 'dn)^] 1 / 2 
the sound velocity of a beam with constant density n\. 
The linearization (8) is valid provided |(5n(x)/ni| 
k\/c\. In this regime, if one further imposes v\ » ci, the 
up-stream density oscillations can be analyzed in term of 
incident and reflected particles (and not quasi-particles). 
This allows to unambiguously define the incident, re- 
flected, and transmitted current as Ji = (ni + 8n\/2)v\, 
J r = 5niV\/2 and J t — n\V\ = (where Sn 1 = 
A[E c i — W{ni)\l n\). Hence, once E d is known, the 
transmission at given incident current Ji is determined 
through 
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Ed = W(m) + ^Sm = W(m) + g-J 4 (l - T) . (9) 

The linearization procedure explained so far is valid in 
the case of small upstream interaction (this is the essence 
of the condition v\ Ci). However, all interaction effects 
are fully taken into account in the down-stream region, 
where they are indeed more important (the constriction 
acts as a barrier which lowers the velocity of the down- 
stream flow and, by flux conservation, increases its den- 
sity [10]). In the following we solve the exact non-linear 
equation (3) and use the linearization procedure only to 
define the transmission coefficient T. Thus, the results 
presented below are of very general validity, but their 
analysis in term of transmission coefficient is only cor- 
rect in so far as the linearization procedure is valid. 

The method is now the following: for a given /x and Jj, 
assume a particular value of T. This determines Joe = 
T Ji, fixes the form of the function Wo(n), the value of 
n(+oo) (it is equal to m,o in case A and to 77,2,0 in case 
B) and of E d = Wb[n(+oo)]. Integrating (3) backwards 
from x = +00 to x = —00 yields E c t, which should be 
compatible with (9). If not, the value of T has to be 
modified until self-consistency is achieved. 

To understand the physical picture we consider an 
abrupt step-like constriction. In this geometry, numerical 
integration of (3) can be bypassed because E c i is simply 
expressed in terms of E d (see Eq. (10)). However, the 
adiabatic approximation (1) is based on the assumption 
that the typical longitudinal length scale is much larger 
than the transverse one. This is clearly violated by an 
abrupt constriction. Nevertheless, in certain parameter 
ranges, the adiabatic approximation remains valid. In 
order to illustrate this point we compute the transmis- 
sion of non-interacting atoms [11], for which the exact 
solution can be obtained numerically. We thus consider a 
linear wave moving in a guide with harmonic confinement 
whose transverse pulsation changes abruptly (at x = 
say) from (up-stream) to lu^ = au;^ (down-stream), 
with a > 1. The incident atoms occupy the transverse 
ground state of the up-stream potential. The numerical 
solution of the problem can be worked out by a straight- 
forward 3D generalization of the procedure devised in 
Ref. [12] for studying a similar 2D problem. We denote 
the transmission T L , the superscript recalling that we are 
in a linear (i.e., non-interacting) regime. The result for 
T L is presented in Fig. 2 for a — 2 and a = 3. The 
vertical bars indicate the location of the energies of the 
reflected (thin lines) and transmitted (thick lines) chan- 
nels. Conservation of angular momentum along the lon- 
gitudinal axis imposes selection rules between channels. 
These rules effectively forbid half of the energetically al- 
lowed channels, and those henceforth do not play any role 
in the transmission. 

Within the adiabatic approximation, the constriction 
is described by an abrupt longitudinal step potential of 



height Vo = lo^ — lu^ . The corresponding transmission 

is Tadia = K,)] 1/2 (VM + V^ - Vo)- 2 (represented 

by a dashed curve in Fig. 2). The onset of transmission 
occurs at fi = Vo, i.e., /j,/u>± = ct — 1. One notices in 
the figure that when /i is increased from this value, the 
adiabatic approximation is initially quite accurate. How- 
ever, at larger values of /i deviations from adiabaticity are 
clearly visible. For a = 2, a sudden lowering of T L oc- 
curs when a new reflected channel opens at /i/a/J; = 2 (in 
all the following we denote as "open channels" those al- 
lowed by energy conservation and symmetry rules). From 
there on, T L diminishes until a new transmission channel 
opens, at /i/w^ = 5. At this point, a sudden increase of 
T L is observed. At large values of /i, T L tends to unity, 
as it should. For a = 3, the process is similar, but the 
breakdown of the adiabatic approximation occurs earlier 
because the opening of the initial transmitted channel - 
at fi/uj^ = 2 - coincides with that of the first allowed 
excited reflected channel. 

In the following we will concentrate on the region 
li/u)^ > a — 1 where the adiabatic approximation is well 
justified and where, as we shall now see, non-linear effects 
may induce strong modifications of the transmission. 
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FIG. 2. Solid curve: transmission T L of a linear wave 
as a function of fj, for a = 2 (top) and a — 3 (bottom). 
Dashed curve: adiabatic approximation. In each plot the ver- 
tical thick (thin) lines indicate the location of the transverse 
eigen- energies in the down-stream (up-stream) confining po- 
tential. 

We now turn to the non-linear problem and consider 
the case V\\{x < 0) = 0, e(n;x < 0) = e(n) and V\\(x > 
0) = Vo, e(n;x > 0) = eo(n). Eq. (3) admits the first 
integral (4) for all x < and (5) for all x > 0. E c i is 
determined through E d by imposing continuity of A and 
A' at x = 0: 

E cl - W[n(0)} - E° cl - W [n(0)} . (10) 

Let's consider case A first. The asymptotic down-stream 
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density is rti.o and thus one has, for all x > 0, n(x) = ni ; o 
(the fictitious classical particle remains at the bottom of 
the potential well Wo)- In particular, n(0) = ni ; o and 
the matching (10) determines E c i uniquely. The value of 
T is denoted T A in this case. 

Case B is more interesting because the structure of the 
down-stream solution is richer: n(x > 0) being part of 
the profile of a gray soliton, n(0) can be varied contin- 
uously provided the matching (10) is fulfilled at a value 
acceptable for Eq. (4). Effectively, the only restriction 
imposed is that n min < n(0) < n max (n min and n max 
are defined in Fig. 1). As a result, for fixed \x and Jj, E c i 
is not uniquely determined by n(+oo), and the transmis- 
sion T B varies between and a value that we denote as 

1 max' 

To be specific, we consider a continuous beam of 23 Na 
atoms propagating through a guide with a transverse con- 
finement u< = 2tt x 2 kHz (in the region x < 0), to which 
we impose a narrowing u>j_ = 2ir x 6 kHz (in the region 
x > 0). This represents a barrier of height Vb = 192 
nK. In the non-interacting case the transmission T L as 
a function of \i is plotted in the bottom part of Fig. 2 
(a = 3). We take \x = 210 nK (this corresponds to the 
kinetic energy of atoms having a velocity of 1.2 cm/s). 
This value of \x corresponds, in the non-interacting case, 
to a regime where the adiabatic approximation holds and 
yields a transmission T L ~ T^ dia ~ 0.70. The repulsive 
interaction between atoms introduces in Eq. (3) a non- 
linear term which, in the region x < 0, reads e(n) = gn 
with g — 2a sc u>< ~ 530 nK.nm. In the region x > 0, 
the transverse frequency of the guide is multiplied by 3, 
and thus eo(n) — 9o n with go = 3 x g. An important 
parameter of the system is the maximum transmitted 
current J^ ax above which no stationary flow can exist in 
the down-stream part of the guide. From Eq. (7) (with 

= 1) one obtains J^ ax — 1.6 x 10 4 atom/s. 

Fig. 3 summarizes the results obtained. The lineariza- 
tion condition v\ 3> c\ is extremely well satisfied in the 
whole range of incident currents considered (the less fa- 
vorable case occurs at large Jj, where v\ ~ 25 c\). The 
horizontal dashed line is the value of the transmission 
coefficient of non- interacting atoms : T L ~ 0.7. It is cur- 
rent independent. When the current is increased from 
zero, T A decreases from this value down to T A ~ 0.5. 
At this point (located with a black spot on the figure), 
Joc(= T A Ji) is equal to J^ ax , and a stationary flow of 
type A is no longer permitted (actually it bifurcates to 
a type B solution). The prominent feature of the behav- 
ior of T A as a function of Jj is its decrease compared 
to the non- interacting value T L . The physical reason 
behind this phenomenon is simple: the available kinetic 
energy necessary to step over the barrier is reduced when 
the interaction energy increases, i.e., when the incident 
current increases. This picture is supported by a pertur- 
bative treatment which accurately describes the flow at 
low incident current (Jj <C J?£ ax ) and confirms that the 



decrease of T A corresponds to an increased fraction of 
the interaction energy in the chemical potential. 




FIG. 3. T A and T B ax as a function of Ji/J™ ax (at fixed 
fi = 210 nK). The horizontal dashed line is the transmission 
T L of non-interacting atoms. The gray zone above the dashed 
hyperbola (of equation T = J^ ax / Ji) is a region where no 
stationary flow exists. Inset: Schematic of the constriction's 
geometry. 

Case B being mediated via interaction, does not ex- 
ist for low current. It exists only above a critical current 
(2g )- 1 (f J ,-V ) 2 (8(j,)- 1 / 2 ^ 1.5 xlO 3 atom/s ~ 0.1 J™ ox . 
From this point, T^ ax increases rapidly up to 1 (reached 
at Jj ~ 0.45 J™ ax in the case of Fig. 3), and then de- 
creases down to a point where one can show that it ex- 
actly meets the end point of T A . From there on, the value 
of T B is limited by the condition that the flow should be 
stationary, and one has T^ ax — J™ ax /Ji, which coin- 
cides with the dashed hyperbola in Fig. 3. We emphasize 
that stationary solutions of type B with arbitrary trans- 
mission < T B < T^ ax exists for any current above the 
critical one. 

The nonlinear transport induced by the repulsive two- 
body interaction has therefore a non-trivial consequence: 
new solutions - of solitonic character - emerge; they allow 
for an increased transmission. This contrasts with the be- 
havior of case A where the transmission is lowered by the 
interaction. One can show that, for any value of /i > Vb, 
there always exists a value of Jj such that complete trans- 
mission exists in case B. The profile for T^ ax — 1 consists 
of a constant up-stream density n(x < 0) = n\ connected 
at x = to half a soliton. Fig. 4 displays the density 
profiles of two stationary flows (case A and case B at 
T B = T r B aax = 1) at an incident current J = 0.36 J£ ax 
(see Fig. 3). Note the significant difference in the densi- 
ties of the up-stream and down-stream profiles, a purely 
nonlinear effect mediated by a solitonic profile. We have 
performed numerical computations that show that the 
same type of solution also exists for smooth constrictions 
and that they are dynamically stable, as confirmed by a 
Bogoliubov analysis. The interactions can thus have two 
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different and, in some sense, opposite consequences on 
the transport properties of a condensate flow. They di- 
minish the transmission in some instances (case A), but 
also allow for new stationary flows that can be perfectly 
transmitted (case B) [13]. 
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FIG. 4. Density profiles (in dimensionless units) for 
the constriction denned in the text, at incident current 
Ji = 0.36 J™ ax (flowing from left to right) and chemical 
potential /i = 210 nK. Thick solid line: solitonic flow with 
T B = T£ ax = 1; thin solid line: type-A flow with T A = 0.7. 

We have restricted our analysis to stationary config- 
urations, but dynamical effects are certainly of interest. 
Amongst these, one could address the question of the 
transient that exists before stationarity is reached, or the 
nature of the flow in parameter regions where stationarity 
is not possible. An other open problem, that clearly de- 
serves investigation, is the dynamical selection of the dif- 
ferent types of stationary profiles discussed above: given 
an initial low density flow (which, from Fig. 3, is of type 
A), which branch (A or B) will be followed when the 
incident flux increases ? 

There are different ways to experimentally realize the 
effect discussed in the present work, namely enhanced 
solitonic transmission of matter waves. We have studied 
one possible implementation, where the step-like poten- 
tial in the longitudinal motion of the condensate is pro- 
duced by a constriction of the guide. An other possibility 
is to apply a blue-detuned laser beam on the region x > 
of a condensate propagating along a guide of constant di- 
ameter. In this case, the characteristics of the flow and 
the barrier should be easily controlled by modifying the 
laser's frequency, intensity and waist, thus allowing for a 
neater experimental observation of the above predicted 
transport phenomena [14]. 
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